Search for: All records

The physical fidelity of turbulence models can benefit from a partial resolution of fluctuations, but doing so often comes with an increase in computational cost. To explore this trade-off in the context of wall-bounded flows, this paper introduces a framework for turbulence-resolving integral simulations (TRIS) with the goal of efficiently resolving the largest motions using a two-dimensional, three-component representation of the flow defined by instantaneous wall-normal integrals of velocity and pressure. Self-sustaining turbulence with qualitatively realistic large-scale structures is demonstrated for TRIS on an open-channel (half-channel) flow configuration using moment-of-momentum integral equations derived from Navier–Stokes with relatively simple closure approximations. Evidence from direct numerical simulations (DNS) suggests that TRIS can theoretically resolve$$35\,\%{-}40\,\%$$of the turbulent skin friction enhancement for friction Reynolds numbers between$$180$$and$$5200$$, without a noticeable decrease or increase as a function of Reynolds number. The current implementation of TRIS can match this resolution while simulating one flow through time in$${\sim}1$$minute on a single processor, even for very large Reynolds numbers. The framework facilitates a detailed apples-to-apples comparison of predicted statistics against data from DNS. Comparisons at friction Reynolds numbers of$$395$$and$$590$$show that TRIS generates a relatively accurate representation of the flow, while highlighting discrepancies that demonstrate a need for improving the closure models. The present results for open-channel flow represent a proof of concept for TRIS as a new approach for wall-bounded turbulence modelling, motivating extension to more general flow configurations such as boundary layers on immersed objects. 

Turbulent boundary layers on immersed objects can be significantly altered by the pressure gradients imposed by the flow outside the boundary layer. The interaction of turbulence and pressure gradients can lead to complex phenomena such as relaminarization, history effects and flow separation. The angular momentum integral (AMI) equation (Elnahhas & Johnson,J. Fluid Mech., vol. 940, 2022, A36) is extended and applied to high-fidelity simulation datasets of non-zero pressure gradient turbulent boundary layers. The AMI equation provides an exact mathematical equation for quantifying how turbulence, free-stream pressure gradients and other effects alter the skin friction coefficient relative to a baseline laminar boundary layer solution. The datasets explored include flat-plate boundary layers with nearly constant adverse pressure gradients, a boundary layer over the suction surface of a two-dimensional NACA 4412 airfoil and flow over a two-dimensional Gaussian bump. Application of the AMI equation to these datasets maps out the similarities and differences in how boundary layers interact with favourable and adverse pressure gradients in various scenarios. Further, the fractional contribution of the pressure gradient to skin friction attenuation in adverse-pressure-gradient boundary layers appears in the AMI equation as a new Clauser-like parameter with some advantages for understanding similarities and differences related to upstream history effects. The results highlight the applicability of the integral-based analysis to provide quantitative, interpretable assessments of complex boundary layer physics. 

Understanding and predicting turbulent flow phenomena remain a challenge for both theory and applications. The nonlinear and nonlocal character of small-scale turbulence can be comprehensively described in terms of the velocity gradients, which determine fundamental quantities like dissipation, enstrophy, and the small-scale topology of turbulence. The dynamical equation for the velocity gradient succinctly encapsulates the nonlinear physics of turbulence; it offers an intuitive description of a host of turbulence phenomena and enables establishing connections between turbulent dynamics, statistics, and flow structure. The consideration of filtered velocity gradients enriches this view to express the multiscale aspects of nonlinearity and flow structure in a formulation directly applicable to large-eddy simulations. Driven by theoretical advances together with growing computational and experimental capabilities, recent activities in this area have elucidated key aspects of turbulence physics and advanced modeling capabilities. 

The data is from a direct numerical simulation of forced isotropic turbulence on a 10243 periodic grid, using a pseudo-spectral parallel code. Time integration of the viscous term is done analytically using integrating factor. The other terms are integrated using a second-order Adams-Bashforth scheme and the nonlinear term is written in vorticity form1. The simulation is de-aliased using phase-shift and a 2√2 /3 truncation2,3. Energy is injected by keeping constant the total energy in modes such that their wave-number magnitude is less or equal to 2. After the simulation has reached a statistical stationary state, 5028 frames of data, which includes the 3 components of the velocity vector and the pressure, are generated and ingested into the database. The duration of the stored data is about five large-eddy turnover times.